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Decaying dark matter particles could be indirectly detected as an excess over a simple power law 
in the energy spectrum of the diffuse extragalactic gamma-ray background. Furthermore, since the 
Earth is not located at the center of the Galactic dark matter halo, the exotic contribution from dark 
matter decay to the diffuse gamma-ray flux is expected to be anisotropic, offering a complementary 
method for the indirect search for decaying dark matter particles. In this paper we discuss in detail 
the expected dipole-like anisotropies in the dark matter signal, taking also into account the radiation 
from inverse Compton scattering of electrons and positrons from dark matter decay. A different 
source for anisotropies in the gamma-ray flux are the dark matter density fluctuations on cosmic 
scales. We calculate the corresponding angular power spectrum of the gamma-ray flux and comment 
on observational prospects. Finally, we calculate the expected anisotropies for the decaying dark 
matter scenarios that can reproduce the electron/positron excesses reported by PAMELA and the 
Fermi LAT, and we estimate the prospects for detecting the predicted gamma-ray anisotropy in the 
near future. 



I. INTRODUCTION 

The full-sky observations of gamma rays undertaken 
in the 1990s by the Energetic Gamma Ray Experiment 
Telescope (EGRET) revealed a map which is highly 
anisotropic. It showed a number of resolved sources [l|, 
such as blazars, as well as an unresolved component, 
which can be attributed almost entirely to Galactic emis- 
sion. The mechanisms which produce the diffuse Galac- 
tic emission (inverse Compton scattering, tt*^ production 
and bremsstrahlung) are well understood. However, the 
actual intensity of each of these contributions depends 
on many details of the propagation of cosmic rays in the 
Galaxy, including the interstellar radiation and magnetic 
fields and the Galactic gas distribution. On the other 
hand, over the last years a picture of cosmic-ray propa- 
gation has emerged which can account for the observed 
abundances of almost all secondary cosmic-ray species. 
Remarkably, this same propagation model can reproduce, 
using the interstellar gas and radiation field distributions 
and a Galactic magnetic field inferred from observations, 
the full-sky gamma-ray map with rather good accuracy. 
However, this requires the introduction of an additional, 
a priori undetermined, diffuse component which is postu- 
lated to be of extragalactic origin and thus isotropic. Its 
energy spectrum is determined by observations at high 
Galactic latitudes. 

Recently it has become apparent that the state-of-the- 
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art propagation models fail to reproduce the measure- 
ments of the positron fraction at energies larger than 
10 GeV. Namely, the secondary positron flux calculated 
with current propagation models, together with the total 
electron-plus-positron spectrum measured by the Fermi 
LAT, yields a positron fraction which monotonically de- 
creases with energy, whereas the PAMELA observations 
reveal that above 10 GeV the positron fraction increases 
with energy @, [3, 0]- This puzzling behavior, already 
suggested in the past by a series of experiments such as 
HEAT [F] , has been interpreted as evidence for a primary 
component of the positron flux (possibly accompanied 
by an identical electron flux). Furthermore, the conven- 
tional propagation model adopted by the Fermi collabo- 
ration, dubbed "model 0" , which fits well the low-energy 
data points of the total electron-plus-positron flux and 
the positron fraction, fails to reproduce the observations 
at higher energies, reinforcing the necessity of an exotic 
source of primary electrons and positrons. 

The most common astrophysical explanation of the 
electron/positron excesses is the electron-positron pair 
production by the interactions of high-energy photons 
in the strong magnetic field of nearby pulsars, such as 
Geminga or Monogem 0, 0, H, [^. However, this inter- 
pretation requires a rather large percentage of the total 
spin-down power injected in the form of electron-positron 
pairs, namely about 40%, and a large cutoff of the elec- 
tron/positron energy spectrum at about 1 TeV. Alterna- 
tively, the electron/positron excesses could be explained 
by the combined emission of both nearby and distant 
pulsars, this solution requiring a percentage of spin-down 
power ranging between 10 — 30% and again a large cutoff 
in the energy spectrum, 800 - 1400 GeV 0]. 

An arguably more exciting explanation of the cosmic- 
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ray electron/positron excesses is the possibility that the 
electrons and positrons are produced in the annihila- 
tion O, [nl, M, [ll or the decay [13, M, M, [13, [11 of 
dark matter particles. Despite the overwhelming gravi- 
tational evidence for the existence of non-baryonic dark 
matter, which makes up more than one- fifth of the energy 
budget of the Universe [l9l| . v erv little is known about its 
particle nature (see Ref. [20| for a review). The most 
popular type of dark matter candidate, the weakly inter- 
acting massive particle (WIMP), can naturally reproduce 
the observed dark matter abundance due to effective self- 
annihilation in the early Universe, after being in thermal 
equilibrium with the baryons before. Today, this same 
annihilation process is expected to produce a possibly ob- 
servable contribution to the measured cosmic-ray fluxes 
on Earth. Detection of such an indirect signal would 
be the first non-gravitational evidence for dark matter, 
with paramount importance to the understanding of its 
nature. 

A lot of effort has been made to study the prospects 
and predictions of cosmic-ray signatures from annihilat- 
ing dark matter, see e.g. [U, [22, 23). However, this is 
not the only possibility for the indirect detection of dark 
matter. Many dark matter models predict that the dark 
matter particle is unstable and decays with very long life- 
times [M [13, [II [13, [11. If the decays occur at a sufli- 
ciently large rate, the decay products could be observable 
as an exotic contribution to the high energy cosmic ray 
fluxes of gamma rays, electrons, positrons, antiprotons, 
neutrinos or antideuterons [H, [13, [H, [H, [10, [Ml [11 . 
Among these, the gamma-ray channel is probably the 
most important to study, due to its sensitivity to far- 
distant sources and its potential to discriminate between 
signals from annihilating or decaying dark matter and 
astrophysical sources. 

In this paper we will present a detailed study of the 
peculiar predictions for gamma rays from decaying dark 
matter. We will concentrate on their angular anisotropics 
on large and small scales. These can be used to discrim- 
inate this component from other contributions to the ex- 
tragalactic diffuse emission. It turns out that the pre- 
dictions for decaying dark matter are much more robust 
than the ones for annihilation, which makes this scenario 
very predictive and easier to confirm or falsify [11 . Sec- 
ondly, we will discuss the prospects to see such signals in 
the upcoming Fermi LAT gamma-ray data, and we will 
apply our results to the decaying dark matter explana- 
tion of the electron/positron excess. 

This paper is organized as follows: In section [TT] we 
will review the basic concepts about production, propa- 
gation and absorption of gamma rays from dark matter 
decay. In section IIIII we will calculate the gamma-ray 
anisotropy expected from the decay of dark matter par- 
ticles on large angular scales. We will show in particular 
the expected anisotropy in scenarios which can explain 
the electron/positron excesses observed by PAMELA and 
the Fermi LAT, and we will argue that this anisotropy 
should be seen by the Fermi LAT. In section [TVl we will 



calculate the angular power spectrum of gamma rays 
from decaying dark matter on small angular scales and 
comment on its observational prospects. Lastly, in sec- 
tion |Vl we will present our conclusions. We also present 
two appendices in which wc discuss statistical proper- 
ties of the large-scale anisotropy and the general obser- 
vational strategy for gamma rays from dark matter decay. 



II. GAMMA RAYS FROM DARK MATTER 
DECAY 

The total gamma-ray flux from dark matter decay re- 
ceives several contributions. The first one stems from the 
decay of dark matter particles in the Milky Way halo and 
reads 



dJh 



alo 



1 



47r mdm Tdm dE, 



7 -'0 



ds Phaio[r{sJ,b)] , 
(1) 

where dN^/dE^ is the energy spectrum of gamma rays 
produced in the decay of a dark matter particle and 
Phaio('") is the density profile of dark matter particles in 
our Galaxy, as a function of the distance from the Galac- 
tic center, r. The received gamma-ray flux depends on 
the Galactic coordinates, longitude I and latitude 6, and 
is given by a line-of-sight integral over the parameter s, 
which is related to r by 



2sR(:) COS b cos / 



(2) 



Here, Rq — 8.3 kpc denotes the distance of the Sun to 
the Galactic center (34] . 

The latest A^-body simulations favor the Einasto den- 
sity profile m, [H, [H 



P^:ir°{r) « exp 



(3) 



which we use throughout the work when not stated oth- 
erwise, and for which we adopt a = 0.17 and the scale 
radius = 20 kpc. For comparison we will also show 
results for the much shallower isothermal profile 
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(4) 



with Ts = 3.5 kpc. We use the local dark matter density 
as determined in Ref. [H, = 0.385 GeVcm"'^, to 
normalize the profiles at the position of the Sun. Related 
uncertainties and their impact on our results will be 
discussed below. 

The electrons and positrons that can be produced in 
the decay of dark matter particles, and which may be 
the origin of the PAMELA and Fermi LAT anomalies, 
also generate a contribution to the total gamma-ray flux 
through their inverse Compton scattering on the inter- 
stellar radiation field (ISRF), which includes the CMB, 
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thermal dust radiation and starlight. Recently, ICS ra- 
diation in connection with the PAMELA excess was dis- 
cussed in Refs. [H, IH, \M, \MJM, El; a pedagogical 
review can be found in Ref. 44]. Furthermore, the in- 
teractions of energetic electrons and positrons with the 
Galactic magnetic field produce synchrotron radiation in 
the radio band with frequencies 0(0.1 — 100 GHz), which 
could also be observed (see e.g. Ref. [Zsl. [46t). 

The production rate of gamma rays with energy 
at the position f of the Galaxy, due to inverse Compton 
scattering of dark matter electrons (or positrons) with 
number density /g± {Ef, , r) on photons of the ISRF with 
number density /iSRF(e,?0; given by 



di?J,^(r) 



dE^ 



de 



' da'^{E,,e) ^ 
dEe — /e±(-Ee,r)/iSRF(e,r) 



(5) 



Here, da^'~^/dE^ denotes the differential cross section of 
inverse Compton scattering of an electron with energy 
Eg, where an ISRF photon with energy e is up-scattered 
to energies between E.y and Ej + dE^ . It can be derived 
from the Klein-Nishina formula and is given by 



da"-{Ee,e) _ 3 ctt 



dE-, 



4 7? 



2q\nq+l + q-2q'^ 



1 (qr)2 



21 + qT 



(1-9) 



(6) 



where ctt — 0.67 barn denotes the Compton scattering 
cross section in the Thomson limit, 7e = Ee/rrie is the 
Lorentz factor of the electron, me = 511 keV is the 
electron mass, and we have defined F = 47ee/me and 
q = E^/T{Ee — Ej). Eq. ([6]) holds in the hmit where 
e,me <C Ee, and kinematics and the neglect of down- 
scattering require that e < E^ < {1/Ee + 1/47^6)^^ = 
E"^'^^. In the calculations we will further assume that 
the photon and electron fields are isotropic; taking into 
account the anisotropy of the photons, which are mainly 
produced in the Galactic disk would give 0(10% — 20%) 
corrections to the ICS fluxes [47|. 

The gamma-ray flux from ICS that is received at Earth 
reads 



halo-IC 

'dE~ 



1 f°° 

{l,b)^2-— ds 



dRF[risJ,b)] 



dE^ 



(7) 



where the factor of 2 takes into account the fact that both 
dark matter electrons and positrons contribute equally to 
the total flux of gamma rays. 

For the number density of ISRF photons we will use 
results from Ref. [i^. On the other hand, the number 
density of electrons and positrons from dark matter decay 
follows from solving the appropriate transport equation, 
which incorporates the effects of diffusion, reacceleration 
and convection in the Galactic magnetic field and energy 



losses due to synchrotron emission and inverse Comp- 
ton scattering on the ISRF (see e.g. (ii, [13, [Ml ) . How- 
ever, at higher energies above a few 10 GeV the transport 
equation is dominated by the energy loss terms, and the 
number density of electrons and positrons can be approx- 
imated by 



/e±(i?e,rO 



1 



^1^10(7^) 



'^dm Tdm J E, 



dEe — ^ ■ 8 

dEe ^ ^ 



Here, b{Ec, r) accounts for the energy losses and contains 
a part that comes from ICS on the ISRF, and a part that 
comes from synchrotron losses in the Galactic magnetic 
field, b = 5ics + &syn- We set fe± = outside of the dif- 
fusion zone, which we model by a cylinder of half-height 
L — 3 kpc and radius i? = 20 kpc.^ The impact of a vari- 
ation of the height of the diffusive halo on our results will 
be discussed below. 

The part of the energy loss that is due to ICS is given 

by 



^ICs(£^e,^) 



de 



_ _ , da'^{Ee,e) ^ (9) 
dE^ [E^ - e) — /isRF(e, r) . 



For electron energies Ee = I GeV, 6ics ranges between 
4.1 X lO^i'^GeVs-i and 1.9 x 10-i^GeVs-\ depending 
on r. At higher energies 6ics approximately scales like 
~ E^. On the other hand, the synchrotron loss part reads 



bsyn{Ee,r) = -cttt; 



,52 



(10) 



where B'^/2 is the energy density of the Galactic mag- 
netic field, and we set B = 6/xGexp(— |z|/5kpc — 
r/20kpc) for definiteness [i^. At position of the 
Sun this yields an energy loss of bgyn — 4.0 x 
10-i^(£;e/GeV)2GeVs-i. 

Substituting Eq. ([5]) into Eq. (O, the flux from inverse 
Compton scattering of electrons and positrons from dark 
matter on the Galactic radiation field can be calculated. 

In addition to the gamma-ray fluxes that originate 
from the decay of dark matter particles in the Milky Way 
halo, there exists a largely isotropic contribution gener- 
ated by the decay of dark matter particles at cosmologi- 
cal distances. Analogously to the Milky Way component, 
the latter receives contributions from the direct decay of 
dark matter particles into photons, and from the gamma 
rays produced by the inverse Compton scattering of dark 



^ For some sample decay channels we have cross-checked with Gal- 
PROP v50p, using appropriately modified versions of the model 
50p_599278 (which adopts a diffusive halo with L = 4 kpc) and 
of the annihilation package, that our approximations give correct 
ICS gamma-ray sky maps at the 30% level everywhere in the sky 
for gamma-ray energies above 1 GeV. 
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matter electrons and positrons on the intergalactic radi- 
ation field. 

The direct decay of dark matter particles at cosmolog- 
ical distances produces a gamma-ray flux that is given 
by 

— — - — X 

in'^dm \ 



where H{z) = i?o\/^A + ^m{z + 1)'"^ is the Hubble ex- 
pansion rate as a function of redshift z, and pc — 5.5 x 
10~^ GeV/ cm*^ denotes the critical density of the Uni- 
verse. Throughout this work we assume a ACDM cosmol- 
ogy with parameters JIa = 0.74, fim = 0.26, Sldm = 0.21 
and h = 7?o/100km Mpc~^ = 0.72, as derived from 
the five-year WMAP data tH. 

The extragalactic ICS radiation of electrons and 
positrons from dark matter is expected to come mainly 
from scattering on the CMB [33|. Note that there is 
a similar component that comes from electrons and 
positrons produced in the Milky Way halo, but outside 
of the diffusion zone (see Ref. [i^ for the analogous 
case of annihilation). For dark matter particles with 
masses below 3 — 5 TeV these components are all ex- 
pected to become relevant only for gamma-ray energies 
Ej < 10 GeV. We will neglected extragalactic ICS radi- 
ation from dark matter decay in the present work. Note, 
however, that this radiation would be largely isotropic 
and would somewhat reduce the overall anisotropy of 
the ICS radiation at these lower energies. 

In Eq. pTjl we included an attenuation factor for 
the gamma-ray fiux, which incorporates the effects of 
electron-positron pair production by collisions of gamma 
rays from dark matter decay with the extragalactic back- 
ground light emitted by galaxies in the ultraviolet, optical 
and infrared frequencies [s^ . The attenuation factor is 
determined by the optical depth t(Ej,z), for which we 
will use the results from [SS^ throughout this work.^ In 
Fig. [T]we show isocontours of the optical depth in the red- 
shift vs. energy plane. It is apparent from the plot that 
gamma rays with high energies around 1 TeV are strongly 
attenuated and come mainly from redshifts z < 0.05. On 
the other hand, the flux of gamma rays originating from 
the decay of dark matter particles in the Galactic halo 
is barely attenuated by pair production on the ISRF at 
energies below 10 TeV (55]. 

To compare the sizes of the halo and the extragalac- 
tic component of gamma rays from decaying dark mat- 
ter (neglecting ICS radiation for simplicity), we show in 



^ In Ref. [53|| the optical depth is calculated for redshifts z < 5. 
Following Ref. [5J|i we assume that the optical depth does not 
increase beyond 2 = 5 and set T{E-y)zy5 = T(E-y)z—5. 
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FIG. 1: Isocontours of the optical depth t{E^,z) of gamma- 
ray photons, emitted at redshift z and observed at Earth with 
energy E-y. We show results for two different models of the 
intergalactic background light [H^l ■ Throughout this work, we 
will adopt the "fast evolution" model. 



Fig. [5] the total flux of photons, integrated over all en- 
ergies, as a function of the angular distance ip from the 
Galactic center, for a dark matter particle which decays 
producing a monoenergetic photon with an energy in the 
range E!y ~ 10 GeV— 1 TeV. As apparent from the figure, 
the cosmological contributions decrease with energy due 
to the attenuation described above, while the radiation 
profile from decaying particles in the halo is independent 
of energy. The halo contribution typically dominates the 
total flux independently of the halo profile, except at low 
energies E!^ < 10 GeV in the direction of the Galactic 
anticenter. The differences in the two dark matter pro- 
files become only relevant near the Galactic center when 
"0 ^ 10°, and at the Galactic center the fiux predicted 
for the Einasto profile is almost one order of magnitude 
larger than the corresponding fiux from the isothermal 
profile. 



III. LARGE-SCALE ANISOTROPIES IN 
GAMMA RAYS FROM DECAYING DARK 
MATTER 

The decay of dark matter particles can produce gamma 
rays that could be detected as an exotic contribu- 
tion to the diffuse extragalactic gamma-ray background 
(EGBG). The diffuse extragalactic background at high 
energies is believed to be dominated by the emission from 
unresolved active galactic nuclei and is expected to ap- 
proximately follow a simple power law, with an intensity 
and index that has to be determined by fitting to the 
data (5^,^3,111]. Thus, if dark matter particles decay at 
a sufficiently fast rate, one generically expects to observe 
a deviation from a simple power law in the gamma-ray 
energy spectrum, which could show up in experiments 
like Fermi LAT. 
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FIG. 2: Angular profile of the gamma-ray signal from dark 
matter decay as function of the angle ip to the center of the 
galaxy. The solid [dashed) line shows the contribution from 
decay in the Milky Way halo, assuming the Einasto (isother- 
mal) profile. Extragalactic contributions are shown in dotted 
lines for the case that dark matter decay produces a monoen- 
ergetic line with energies between E'^ = 10 and 1000 GeV. 
The fluxes are integrated over energy and normalized to the 
size of the extragalactic component when absorption is ne- 
glected. 



A complementary signature of dark matter decay 
is the observation of anisotropics in the EGBG. It is 
well known that the offset between Sun and Galactic 
center causes a peculiar angular dependence in the 
gamma-ray signal from dark matter decaying [ssj (or 
annihilating [sl, [b^I) in the Milky Way halo, even if 
the halo profile itself is isotropic. The halo signal is 
largest in direction of the Galactic center and smallest 
in direction of the Galactic anticenter (c/. Fig. [2]). The 
observation of an anisotropy that is aligned in this way 
would be a strong signal for a contribution from dark 
matter, and, on the other hand, its non-observation 
would provide strong constraints. Gamma rays from the 
decay of dark matter particles at cosmological distances 
are isotropic and tend to reduce the anisotropy of the 
overall signal. This attenuation effect is however always 
small, due to the relative weakness of the extragalactic 
component. 

To analyze the prospects of detecting a gamma-ray 
anisotropy from dark matter decay at the Fermi LAT 
it is convenient to define the quantity 



A 



Jgc — Jgac 
Jgc + Jgac 



(12) 



where Jgc and Jgac in general denote the total diffuse 
gamma-ray flux (from dark matter and from astrophysi- 
cal sources) integrated over E-y in some energy range, and 
averaged over the hemisphere in direction of the Galac- 
tic center (GC) and anticenter (GAG), respectively. Sky 
regions with small, |&| < bo, or large, |6| > 6i, Galactic 



Sky patch 


Anisotropy A(,o:i,i 


6o 


6i 


Einasto 


Isothermal 


10° 


90° 


0.21 -0.36 


0.20 - 0.33 


10° 


20° 


0.32 - 0.50 


0.29 - 0.45 


20° 


60° 


0.21 -0.35 


0.20 - 0.33 


60° 


90° 


0.07-0.13 


0.07-0.13 



TABLE I: Anisotropy of the gamma-ray signal from the de- 
cay t/j — > of a fermionic dark matter particle ^p, after sub- 
tracting astrophysical contributions. The ranges correspond 
to the anisotropics from gamma-ray lines with energies be- 
tween 10 GeV and 1000 GeV in different regions of the sky, 
see Eq. 

As an example, we calculate the anisotropy parameter 
A in different regions of the sky for the decay ij) 
of a fermionic dark matter particle ip, after subtracting 
astrophysical sources (the Galactic foreground and 
the extragalactic gamma-ray background), in order 
to compare the anisotropy expected purely from dark 
matter decay to the anisotropy expected from the 
Galactic models (in the rest of the paper, however, we 
will consider both sources of gamma rays simultaneously, 
calculating the anisotropy of the total flux). In this case, 
the energy spectrum of gamma rays has two components: 
a monoenergetic line from the decay of dark matter 
particles in the halo and a redshifted line from decays at 
cosmological distances (note that in this decay channel 
there is no contribution from IGS). We show in Tab. U 
the values of A for dark matter particles with masses 
between 20 GeV and 2 TeV, producing monoenergetic 
photons with energies 10 GeV and 1 TeV, respectively. 
It is interesting that the anisotropy parameter can be as 
large as 0.5 for large energies and relatively low latitudes. 
In the region defined by 6o — 10° and hi — 90° (on 
which we will concentrate below), the anisotropy of the 
"pure" dark matter signal ranges between 0.20 and 0.36, 
with only little dependence on the profile of the dark 
matter halo. These values have to be compared with the 
anisotropies of the Galactic foreground as predicted by 
Galprop (see below), which are considerably smaller, 
and typically A < 0.10 in all the regions that are shown 
in Tab. HI up to energies above 300 GeV. Furthermore, 
the anisotropies measured by EGRET for energies 
below 10 GeV are consistent with the predictions for the 
Galactic foreground [6l|. 

From the theoretical point of view, the search for 
anisotropies in the gamma-ray fiux is a cleaner method 
for the indirect detection of dark matter than the 
search for an excess in the spectrum of the EGBG. As 
mentioned above, the genuinely extragalactic flux from 
active galactic nuclei and other extragalactic sources is 
very poorly understood. Thus, it is difficult to make firm 
predictions for the total gamma-ray flux in scenarios 
with decaying dark matter, even when the particle 



6 



physics model is specified (namely, the dark matter 
mass, lifetime and decay modes). Moreover, there are 
other potentially important isotropic contributions to 
the total flux with an intensity that cannot be predicted 
theoretically. For instance, interactions of high energy 
cosmic rays with debris in the hypothetical Oort cloud 
could produce a sizable gamma-ray flux, provided 
that the column density is larger than 10~'^gcm~^ 
[B^ I- Since all these contributions to the total flux are 
perfectly isotropic, they cancel out when calculating the 
difference of the fluxes between the Galactic center and 
the Galactic anticenter regions. 

To illustrate the large-scale anisotropy A that could be 
produced by dark matter decay, we will show predictions 
for different dark matter decay channels and masses, in- 
cluding Galactic foreground radiation and the astrophys- 
ical extragalactic background. As Galactic foreground we 
take predictions of the conventional GALPROP-model as 
presented in Ref. (model 44_500180), and as astro- 
physical EGBG we will use the following parameteriza- 
tion, which agrees with preliminary results from Fermi 
LAT Q 

(13) 

We show in the left panels of Fig. [3] the predicted 
anisotropy of the total gamma-ray flux that would be 
measurable in different regions of the sky if the dark 
matter particles decay exclusively into r+r^ pairs. The 
anisotropy is calculated with taking into account the 
galactic foreground and the extragalactic gamma-ray 
background as discussed above. The dark matter mass is 
taken to be mdm = 600 GeV, and we set the lifetime to 
Tdm — 3.5 X 10^^ s. The energy spectra dN^^c/dE^^c of the 
photons, electrons and positrons produced in the d ecay 
are calculated with the event generator Pythia 6.4 [6a |. 
The lifetime is chosen such that the gamma-ray fluxes 
are below and compatible with the EGBG, as demon- 
strated in the right panels of Fig. [31 where we also show 
preliminary results from the Fermi LAT for comparison.^ 
Furthermore, we note that the contribution to the local 
electron and positron fluxes from dark matter decay is 
negligible in this scenario. Interestingly, for the adopted 
choice of parameters, an anisotropy is predicted that is 
significantly different to the one expected from the diffuse 
Galactic emission in the conventional Galprop model. 
Such an anisotropy should, moreover, be observable by 
Fermi LAT, as illustrated by the boxes in the figure, 
which correspond to our estimates of the one-year and 



five-year statistical errors of Fermi LAT, assuming expo- 
sures of £ = 3x10^° cm^ s |64] and e = 2 x 10^^ cm^ s [SJ] , 
respectively (for a discussion about our calculation of 
the statistical errors see below and appx. As ex- 

pected, the size of the anisotropy is largest at low lat- 
itudes b < 20°, and decreases slowly when considering 
higher latitudes. On the other hand, the statistical error 
is smallest in the sky-patch shown in the upper panels of 
Fig. [3l where the whole sky with |6| > 10° is included. 
The effects of ICS radiation are negligible in the present 
case of decay into r+r" pairs, since the electrons and 
positrons produced in the subsequent decay of the taus 
have only relatively small energies. 

To illustrate the impact of ICS radiation on the 
anisotropy parameter A, we show in Fig.[4]the anisotropy 
of the gamma-ray flux assuming that the dark matter 
particle decays into e~^e~ pairs (with TOdm = 1000 GeV 
and Tdm = 2 X 10^^ s). In this case the dominant source of 
gamma rays is inverse Compton scattering (note that we 
assumed that the decaying dark matter particle has spin 
1; for scalar dark matter particles, helicity suppression 
leads to an enhanced production of final-state radiation 
[6^ . weakening the relative contribution from ICS). For 
reference, we also show the anisotropy that would be 
measurable if ICS radiation were absent (dashed lines in 
left panels of Fig. |4]). Note that in this scenario electrons 
and positrons produced in the dark matter decay give 
a sizeable contribution to the local cosmic-ray fluxes, 
without being in conflict with the PAMELA and Fermi 
LAT data. Again, we find that a sizable anisotropy is 
expected in several patches of the sky. In this case, 
however, the gamma rays relevant for our predictions 
are mainly produced close to the Galactic center, above 
and below the Galactic disk. Hence, the anisotropics 
are relatively weak at higher latitudes |6| > 20°. A 
more detailed discussion about the differences between 
prompt and ICS gamma rays from dark matter decay 
can be found in appx. [Bj 

If dark matter decay is the origin of the excess in the 
positron fraction observed by PAMELA and in the total 
electron-plus-positron flux observed by the Fermi LAT, 
the predicted anisotropics in the gamma-ray flux can 
be quite large. In Fig. O we show our results for the 
anisotropy. A, which is expected to be observed by the 
Fermi LAT if the dark matter particle decays via one 
of the five different decay channels that were found in 
Ref. [IHl to fit well the positron and electron data.^ To 
minimize the statistical errors, we concentrate on the re- 
gion defined by bo = 10° and 6i = 90° (see appx. iBt . 
As apparent from the plots, the predictions for some 
decay channels (namely the decay into W^n^) are al- 



^ Fitting the prehminary Fermi LAT results with a Galactic fore- 
ground model is well beyond the scope of this paper. Hence, 
there is a mismatch between the total fluxes and the data, which 
does not affect our conclusions. 



* We use slightly larger lifetimes than the ones given in Ref. [Tsl 
to account for the larger value for the local dark matter density 
used in this work. 
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FIG. 3; Left panels: Anisotropy of the gamma-ray flux from dark matter decay into t^t~ pairs as a function of energy. The 
dotted line shows the background anisotropy as expected from the Galactic foreground, while the solid line shows the anisotropy 
of signal + background. We also show the signal + background anisotropy neglecting gamma rays from ICS {dashed, overlapping 
with solid line). The boxes show estimates of the statistical errors for one- year and five-year Fermi LAT observations. Right 
panels: Gamma-ray fluxes averaged over all Galactic longitudes as a function of energy. The thin solid line shows the gamma 
rays from dark matter decay. The two dash-dotted lines show the astrophysical EGBG and the Galactic foreground separately. 
The thick solid line shows the sum of all contributions, whereas the dotted line shows the sum without contributions from 
dark matter. From top to bottom the different panels show predictions for different patches of the sky. The data points show 
preliminary Fermi LAT data 164]] for the total diffuse flux {upper points, in the 10° : 90° region we averaged the presented 
preliminary results appropriately) and the EGBG {lower points). 
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ready in conflict with the preUminary results of the Fermi 
LAT collaboration for the EGBG, whereas other decay 
channels are marginally consistent (e.g. the decay into 
fJ-~^ IJ-~ pairs). However, even for those channels which are 
compatible with the data, sizeable anisotropics, around 
A ~ 0.2 — 0.3, are predicted at energies ~ 100 GeV. 
This is significantly different from the anisotropy ex- 
pected for the astrophysical foreground. As indicated 
by our estimates of the statistical error bars for one- 
year Fermi LAT data taking, this deviation should be 
clearly visible in the upcoming results for the diffuse 
gamma-ray sky. On the other hand, its non-observation 
would set very strong constraints on the decaying dark 
matter interpretation of the positron excess observed by 
PAMEL A/Fermi LAT. 

Uncertainties in the determination of the above 
large-scale anisotropy come from different sources. 
When neglecting ICS radiation, the prediction of an 
anisotropy between 0.2 and 0.3 in the dark matter signal 
at latitudes |6| > 10° is relatively robust. The main 
sources of uncertainty are the profile of the Milky Way 
dark matter halo and its normalization. As discussed 
above, the dependence on the profile is rather weak 
(c/. Tab. m. Only in case of a much lower value for 
the local dark matter density, say = 0.2GeVcm~'^, 
and only for gamma rays with energies E'^ < 10 GeV, 
the anisotropy can become as small as A ~ 0.15. On 
the other hand, the size and anisotropy of the ICS 
radiation from electrons and positrons originating from 
dark matter decay is plagued by many uncertainties like 
the exact height of the diffusion zone, the distribution 
of the ISRF and the size of the Galactic magnetic field. 
In general the ICS radiation, and hence the overall 
anisotropy of the observed flux, becomes stronger if the 
height of the diffusive halo is increased, but a detailed 
study of these uncertainties is beyond the scope of this 
paper. Note, however, that the large scale anisotropies 
predicted for the decay channels shown in Fig. [5] are 
sizeable even if ICS radiation is neglected (dashed lines 
in left panels). 

Finally we will discuss in a more quantitative way 
the main prospects for the Fermi LAT to detect gamma 
rays from dark matter decay through the observations of 
large-scale anisotropies. To this end we will neglect in- 
verse Compton radiation from the electrons and positrons 
produced in the decay of the dark matter particle, and 
we will assume perfect subtraction of the Galactic fore- 
ground. The remaining flux is then expected to be consti- 
tuted by the isotropic EGBG, which is possibly contam- 
inated by anisotropic radiation from dark matter decay. 
For definiteness, we will assume that the remaining flux 
follows the power law in Eq. 

Provided that a fraction fs of the considered gamma- 
ray photons in a given energy range is due to decaying 
dark matter, the measured anisotropy A is given by 



Here, As and Abg denote the anisotropy of the dark 
matter signal, which can be read off from Tab. [H and 
the anisotropy of the astrophysical background, which is 
Abg = in our case, respectively. A possible detection 
at the 3(T-level requires that 



5 (J A 



A^-A 



bg 



(15) 



where aA denotes the standard deviation of the 
anisotropy A. It depends on the total number of mea- 
sured photons, and can be approximated by aA — 



n: 



-1/2 



The adopted approximation for a a is better than 



10% as long as |^| < 0.3 and a a < 0.2. See appx. Elfor 
a short discussion. The photon number is given by 



(16) 



En 
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A = fsAs + (1 - /.)Abg 



(14) 



where e denotes the experimentally given exposure (see 
above) and Slsky is the solid angle of the observed sky, 
which is given by figky — 0.83 • 47r if the Galactic disk 
with |6| < 10° is excluded. Following Eq. the Fermi 
LAT win detect ~ 3.0 x 10'* {N^ ~ 1.1 x 10^) photons 
with energies E^ > 10 GeV {E^ > 100 GeV) after five 
years of data taking. Taking for definiteness As = 0.3, 
this allows in principle a 3cr-detection of a dark matter 
contamination down to fs — 6% {fs — 30%). Note, 
however, that additional statistical noise and systematic 
uncertainties from point source subtraction and the 
determination of the Galactic foreground are neglected 
and can reduce the sensitivity to fs by factors of order 
one. 



IV. ANGULAR POWER SPECTRUM OF 
GAMMA RAYS FROM DECAYING DARK 
MATTER 

We will now discuss the angular power spectrum 
Ci of gamma rays which stem from the decay of 
dark matte r p articles (for the annihilation case, see 
e.g. [6i,[6i,[2a,[7l|,[2l[2i|). On large angular scales the 
angular power spectrum is expected to be completely 
dominated by the dipole-like asymmetry of the halo 
flux, which is - as discussed above - due to the offset 
between the Sun and the Galactic center. However, on 
smaller scales spatial fluctuations of the dark matter 
density can become relevant. These fluctuations are 
related to subhalos of our own Galactic halo, and to the 
large-scale structure of the dark matter distribution in 
the nearby Universe. The abundance and distribution 
of subhalos and their impact on indirect searches for 
dark matter is still a debated question [H, [ill . Recent 
iV-body simulations suggest that inside a radius of 
r — 100 kpc around the Galactic center, the fraction of 
the halo mass that is bound in subhalos can be as small 
as a few percent [s^. Such a small value would lead to 
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FIG. 4: Same as Fig.|3] but for decay into e^e pairs. Here, in the left panels, the difference in the predictions for the anisotropy 
when including ICS radiation (solid line) and neglecting ICS radiation {dashed line) is clearly visible. 
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a strong suppression of subhalo-related signals in case of 
decaying dark matter. In this work we will concentrate 
on the more robust predictions related to the large-scale 
structures of the dark matter distribution, and we will 
assume a smooth Milky Way halo when calculating the 
angular power spectrum. Furthermore, we neglect ICS 
radiation throughout this section. 

The fluctuations of the dark matter gamma-ray flux 
in a certain energy band can be expanded in spherical 
harmonics Yf™ with coefficients aim according to 

SJin) = J{n) - (J) = (J) aimYimin) • (17) 

im 

Here (J) denotes the mean gamma-ray flux, averaged 
over the whole sky, and the coefficients aim can be cal- 
culated from 



= (J)-W darJJ{n)Y;„,{n) 



(18) 



Our goal is to calculate the angular power spectrum 
Ci = (|a£,np), which can be split up into an extragalactic 
and a halo part according to 



Ci = f'cT + (1 - h^yc, 



halo 



(19) 



Here, /eg denotes the fraction of dark matter photons 
with extragalactic origin. Cross-correlations between 
the halo and the extragalactic component vanish in our 



Following Eq. p8p . the angular power spectrum of the 
halo component, C\^^° , can be calculated by 



dV" sin'i/'-P°(cosV') Jhaio(V') 



(20) 



where the denote the associated Legendre functions, 
>/haio(V') is the angular profile of the halo flux as plotted 
in Fig. [21 and the power spectrum has to be normalized 
such that Cj''^'" = 47r. 

For the calculation of the extragalactic component of 
the angular power spectrum, C^^, it is convenient to de- 
fine the window function 



W{z) 



^dmPc 
47rTOdmTdm J E, 



IK 



[{z + 1)E, 



(21) 

for which holds that (Jcg) = dz W (z) / H (z) 

(cf. Eq. (jlip above). The window function mainly de- 
pends on the energy spectrum of gamma rays from the 
dark matter decay dN-y/dE^, and on the optical depth 
t{E^, z). Following the Limber approximation, the angu- 
lar power spectrum of the extragalactic c omp onent can 
be determined by (for details see e.g. Ref. p^ ) 



Cf^{J,^)-^j^'^W{zfp{k=^-,z 



(22) 



Here, P(fc, z) denotes the power spectrum of the dark 
matter distribution as a function of redshift z, and the 
comoving coordinate r is related to z by dr — dz/H{z). 
The dark matter power spectrum at z = can be de- 
rived from A^-body simulations or analytically in the Halo 
model, and we will adopt results from Ref. [54|. The cor- 
responding redshift-dependence follows from the relation 



P{k,z)^P{k),=o-D{zf 



(23) 



where D{z) oc H(z) dz'{l + z')H{z')~^ is the linear- 
growth factor, normalized such that £'(0) — 1 [t^. In 
principle the above set of equations allow a calculation of 
the angular power spectrum of gamma rays for arbitrary 
gamma-ray spectra. 

In Fig. [H] we show the quantity W{z)/H{z) as a 
function of redshift, assuming that dark matter decay 
produces a monoenergetic line with energies between 
E'.y — 10 GeV and 1 TeV. When normalized appropri- 
ately, W{z)/H{z) ■ dz gives the fraction of the received 
extragalactic dark matter photons that were emitted at 
redshifts z . . . z -\- dz. As evident from the figure, the 
region of the Universe that contributes to the observable 
gamma-ray fluxes becomes larger for lower gamma-ray 
energies. At E'^ = 1 TeV, around 90% of the observable 
photons come from comoving distances r < 400 Mpc, 
whereas at lower gamma-ray energies around 10 GeV 
photons are received that come from distances up to 
8 Gpc. Hence, in the latter case fluctuations in the 
extragalactic dark matter signal due to large-scale 
structures, which typically have sizes up to several 
100 Mpc, are efficiently averaged out, which results in a 
correspondingly small value of C^^. For higher gamma- 
ray energies the probed region shrinks, yielding a larger 
C^®. Note, however, that the fraction of gamma rays 
from the decay of dark matter particles at cosmological 
distances, /eg, presents an opposite energy dependence 
and decreases by an order of magnitude, from 0.42 
to 0.03, when increasing the gamma-ray energy from 
E'^ = 10 GeV to 1 TeV. This behavior counteracts the 
growth of in Eq. 

Our results for the total angular power spectrum 
of gamma rays from decaying dark matter are shown 
in Fig. [71 including the halo and the extragalactic 
component, and assuming that the decay of the dark 
matter particle only produces monoenergetic photons 
with energies between E'^ = 10 GeV and 1 TeV. As 
apparent from the plot, the halo component of Ci 
increases slightly with E'^, which is related to a corre- 
sponding decrease of /eg in Eq. p9)l . The variation of 
the extragalactic component with gamma-ray energy E'^ 
is also on the 0(1) level in most of the energy range, 
except at the highest energies E'^ > 1 TeV. In sum- 
mary, the angular power spectrum looks qualitatively 
similar in the whole considered energy regime. The 
halo component dominates at large angles as long as 
I < 100, whereas the extragalactic component takes over 
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FIG. 5: Predictions for the five darlc matter decay channels that where found to fit the positron excess as observed by 
PAMELA/Fermi. The left panels show predictions for the total measurable anisotropy at latitudes |6| > 10° {solid line) and 
estimates of the corresponding 1-year and 5-year statistical error bars of Fermi LAT. The dotted line is the anisotropy predicted 
by Galprop. In the right panels we show predictions for the averaged fluxes in the same sky region. Line coding is the same 
as in Fig. [3] 
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W{z)/H{z), as a function of redshift z, with arbitrary nor- 
malization. We assume that dark matter decay is producing 
a monoenergetic line with energies between E'^ = 10 GeV and 
1 TeV. For comparison, we also show values for the comoving 
coordinate r{z). 
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FIG. 7: Angular power spectrum of gamma-ray signal from 
decaying dark matter. Shown are the contributions from 
the halo component [dashed line) and the contribution from 
the extragalactic radiation [dotted line, based on results from 
[H^l). The dash- dotted lines show for different gamma-ray en- 
ergies our predictions of the photon noise level after five-years 
Fermi LAT observation. 



at smaller angles with t > 100, yielding values for the 
power spectrum around l)Q/27r ~ 0(10^^-10"^). 

To shortly discuss the observational prospects for 
Fermi LAT, we will simply assume that the total EGBG 
(as given by Eq. ([TS]) ) above a certain energy is due to 
dark matter decay. Then, following Ref. [i^, the la de- 
viation of the power spectrum, averaged over i-hvas with 
size A^, is given by 

= '''' 



where Cn = 47r/skyiV*°V(^'')^ is the photon noise, 
and A'^"* denote the number of measured dark matter 
photons and number of total photons, respectively, /sky is 
the observed fraction of the sky and Wg — exp(— cr^^^/2) 
is the window function of a Gaussian point-spread func- 
tion. For the Fermi LAT we take as beam size (Jfc = 0.1°,^ 
and we set = 0.5.^. We only consider latitudes with 
\h\ > 20°, which implies that /sky = 0.66. In this part 
of the sky the ratio between the total gamma-ray fluxes 
and the EGBG flux is 0(3) (see Fig. [5]), and we wiU 
set N^°^ /N^ — 3 for definiteness. The number of signal 
photons follows from Eqs. pT|) and (jl6|) . 

The resulting photon noise (second term of SCg in 
Eq. ([M]) ) for five-years Fermi LAT observation is shown 
in Fig. [7] by the two dash-dotted lines, taking into ac- 
count only gamma rays with energies > 10 GeV or 
> 100 GeV, respectively. As apparent from these lines, 
the small-scale fluctuations are completely hidden under 
photon noise in the energy range we consider, and hence 
they are not accessible by current instruments like Fermi 
LAT. However, note that this is in contrast togamma- 
ray signals from annihilating dark matter l7ll . 17311 or dif- 
ferent extragalactic astrophysical sources |76l . |77[ , where 
the small-scale anisotropics can be orders of magnitude 
larger. This fact might eventually provide a handle on 
distinguishing the different scenarios, but a detailed in- 
vestigation of the prospects is beyond the scope of this 
paper. 



V. CONCLUSIONS 

Dark matter particles could decay into gamma rays 
at a rate which is sufficiently large as to allow their indi- 
rect detection through an excess over the expected power 
law in the energy spectrum of the diffuse extragalactic 
gamma-ray background. In this work we have discussed a 
complementary way of indirectly detecting unstable dark 
matter particles by exploiting the fact that the Earth is 
not located at the center of the dark matter halo. We 
have discussed the relative size of the extragalactic and 
the halo component of the gamma rays from dark matter 
decay, incorporating the attenuation effects from pair- 
production on the intergalactic background light, and we 
have calculated the dipole-like anisotropy between the 
high latitude gamma-ray flux coming from Galactic cen- 
ter and from Galactic anticenter regions for different de- 
cay channels. We have furthermore demonstrated the 
strong impact of gamma-rays from inverse Compton scat- 
tering of electrons and positrons from dark matter decay 
on the anisotropy signal. We have found that if dark 
matter decay is the correct explanation of the excesses in 
the positron fraction and the total electron-plus-positron 
flux reported by the PAMELA and Fermi LAT coUabo- 



^ See |http: //www-glast ■ slac . Stanford. edu| 
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rations, such an anisotropy in the gamma-ray flux should 
be observed by the Fermi LAT (see Fig. . 

Lastly, we have calculated the angular power spectrum 
of the gamma-ray signal from dark matter decay, which 
exhibits imprints from the large scale structure of the 
Universe (see Fig. [7])- We have shown, however, that 
for gamma ray energies above a few GeV these small 
scale effects are too weak to be observed with present 
instruments. 



level as long as (A) < 0.6 and (Ja < 0.2. For small 
enough anisotropics A, however, the standard deviation 
is just given by cr^ = V with good accuracy, and we 
use this approximation throughout the work. 

APPENDIX B: OBSERVATION OF PROMPT 
AND ICS RADIATION FROM DARK MATTER 
DECAY 
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APPENDIX A: STATISTICAL ERRORS OF THE 
LARGE-SCALE ANISOTROPY 

Here, we shortly discuss the statistical errors of the 
large-scale anisotropy as defined in Eq. (I12p . which 
are due to shot noise. Statistical errors of small-scale 
anisotropics are discussed e.g. in Ref. [78*]. 

The measured anisotropy A and the total number of 
measured photons N are related to the number of pho- 
tons measured in direction of the Galactic center, A^i, 
and anticenter, N2, hy A = {Ni - N2)/{Ni + N2) and 
N ~ Ni + N2. The Ni follow a Poisson distribution with 
mean (Ni) and standard deviation a^. — ^ (N^). Con- 
sidering the propagation of uncertainty, it is straightfor- 
ward to derive that the statistical error of the anisotropy 
is given by 



{N) 



(Al) 



which is expected to hold for small enough (A) ~ A and 
large enough (N) ~ N. 

On the other hand, one can derive the exact probabil- 
ity distribution function of the anisotropy A by starting 
with the above Poisson distributions for the Ni, perform- 
ing an appropriate redefinition of the parameters and in- 
tegrating out the total number of measured photons. The 
result is a function of the mean values (A) and {N) and 
can be written in the compact form 



pdf{A) 



(N) 



2{Ni)l{N2 



l + A 



(A2) 

From this equation one can check, for example, that a 
normal distribution with mean (A) and standard devi- 
ation as in Eq. (|Aip gives correct error bars at the 5% 



In this appendix we will briefly discuss the differences 
between observing ICS radiation from positrons and elec- 
trons from dark matter decay and observing the gamma 
rays that come directly from the decay itself (prompt ra- 
diation like internal bremsstrahlung) by means of signal- 
to-noise and signal-to-background ratios. 

Signal-to-noise ratios quantify the significance of a sig- 
nal against statistical noise. The signal-to-noise ratio 
S/N of a dark matter signal with respect to the back- 
ground is given by 



N 



N 



7.S 



7,s 



N 



(Bl) 



7,bg 



where N-y^s and N^^hg denote the number of detected 
signal and background photons, respectively, that are 
observed in a given sky region AQ and energy band 
Ej = Eq . . .El. Signal photons are all photons from 
dark matter decay, background photons are in princi- 
ple all other observed photons, including the astrophys- 
ical part of the EGBG and the Galactic foreground. 
Since the number of detected photons scales like iV^_,; oc 

Si,nd^SEyEdJ,/dE 



I An ""'^ 
tio in the limit AfJ 



^ = Af2 Ji , the signal-to- noise ra- 
is proportional to 



Js 



N 



oc 



(B2) 



where Js and Jbg denote the appropriately averaged and 
integrated signal and background gamma-ray fluxes, re- 
spectively. 

In the upper panel of Fig. [5] we plot the relative signal- 
to-noise ratio (with arbitrary normalization), assuming 
that the background completely dominates the signal, 
Js ^ Jbg, and neglecting ICS and extragalactic radi- 
ation. As background we take the predictions of the 
conventional Galprop model at E^ = 100 GeV (from 
Ref. [11], see above), but the results do not change 
qualitatively for other energies. The anisotropy of the 
dark matter signal as function of I is well recognizable 
in the plot. Furthermore it is apparent that, from the 
perspective of statistical noise, the sky regions that are 
most sensitive to decaying dark matter signals lie close 
above and below the Galactic center, with |Z| < 25° and 
5° < 1^1 < 35°. 

To determine the best observational strategy in light 
of the systematics that are related to the determination 
of the Galactic foreground, it is more convenient to con- 
sider the signal-to-background ratio Js/Jbg (we assume 
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FIG. 8: Upper panel: Relative signal-to-noise ratio of gamma- 
ray signal from dark matter decay as a function of the Galactic 
longitude I and latitude b, normalized to one at the Galac- 
tic center. Lower panel: Relative signal-to-background ratio 
for the same process. Extragalactic and ICS radiation is ne- 
glected, and as background we take the predictions of the 
conventional Galprop model at 100 GeV. 

that systematic uncertainties scale roughly like ~ Jhg)- 
We show the signal-to-background ratio as a function of 
the Galactic coordinates in the lower panel of Fig. [S] 
Again, the large-scale anisotropy of the dark matter flux 
are clearly visible in the plot. Furthermore, it is apparent 
that concerning systematics the best strategy is to avoid 
regions near the Galactic plane and to observe fluxes only 
at higher latitudes, |6| > 20°. However, as a compro- 
mise between statistical and systematic uncertainties we 



choose to consider the whole region 10° < |&| < 90° in 
most of the present paper. 

In contrast to gamma rays that come directly from 
the dark matter decay, the gamma rays that stem from 
inverse Compton scattering of dark matter positrons or 
electrons on the ISRF are mostly coming from the re- 
gion near the Galactic center. This can be seen in Fig. [9l 
where we plot the signal-to-background ratio of the pure 
ICS signal of dark matter decaying into an e^e~ pair 
(with JTidm = 1 TeV) as a function of the Galactic co- 

signal-to-baclcground ratio of ICS radiation 
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FIG. 9: Relative signal-to-background ratio of pure ICS radi- 
ation from dark matter decay at = 10 GeV. We assume 
that mdm = 1 TeV and that the dark matter is decaying into 
e^e~ pairs. As background we take the fluxes from the con- 
ventional Galprop model. 



ordinates. The gamma-ray energy is E-^ — 10 GeV in 
this plot. As background we again use the predictions 
from the conventional Galprop model. From the figure 
it is apparent that the relative size of the signal peaks at 
regions very close below and above the Galactic center, 
with |;| < 20° and 5° < |6| < 30°. This suggests that 
concentrating the observation on these regions is most 
promising for the search for ICS radiation from dark mat- 
ter decay. However, in light of the large underlying un- 
certainties related to the predictions of ICS radiation we 
will neglect these subtleties, and we consider ICS radia- 
tion only in how far it affects the anisotropics and fluxes 
in the sky regions that are most promising for the search 
of gamma rays coming from the dark matter decay itself. 
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